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Tick salivary glands produce complex cocktails of bioactive molecules that facilitate blood 
feeding and pathogen transmission by modulating host hemostasis, pain/itch responses, 
wound healing, and both innate and adaptive immunity. In this study, cutaneous responses 
at Dermacentor andersoni bite-sites were analyzed using Affymetrix mouse genome 
arrays and histopathology at 12, 48, 96 and 120 h post- infestation (hpi) during primary 
infestations and 120 hpi during secondary infestations. The microarray data suggests: (1) 
chemotaxis of neutrophils, monocytes, and other cell types; (2) production and scavenging 
of reactive oxygen species; and, (3) keratin- based wound healing responses. Histological 
analysis supported the microarray findings. At 12 hpi, a mild inflammatory infiltrate was 
present in the dermis, especially concentrated at the junction between dermal connective 
tissue and underlying adipose tissue. A small lesion was located immediately under the 
hypostome and likely represents the feeding "pool." Surprisingly, at 48 hpi, the number 
of inflammatory cells had not increased from 12 hpi, perhaps mirroring the reduction 
in gene expression seen at this time point. The feeding lesion is very well defined, 
and extravasated erythrocytes are readily evident around the hypostome. By 96 hpi, the 
inflammatory infiltrate has increased dramatically and the feeding lesion appears to have 
moved deeper into the dermis. At 120 hpi, most of the changes at 96 hpi are intensified. 
The infiltrate is very dense, the epidermis is markedly thickened, the feeding lesion is 
poorly defined and the dermal tissue near the hypostome appears to be loosing its normal 
architecture. In conclusion, during D. andersoni feeding infiltration of inflammatory cells 
increases across time concurrent with significant changes in the epidermal and dermal 
compartments near the feeding tick. The importance of changes in the epidermal layer 
in the host response to ticks is not known, however, it is possible the host attempts to 
"slough off" the tick by greatly increasing epithelial cell replication. 
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INTRODUCTION 

Hard (Ixodid) ticks are important vectors of disease world- wide, 
making them a significant threat to public health. The process of 
tick feeding is an extended interplay between the tick, tick salivary 
molecules, and the host response. Tick saliva is a complex mix- 
ture of molecules that has been shown to inhibit a broad range of 
host responses and facilitate pathogen transmission, while host 
responses can, in some cases, inhibit the feeding process and 
block pathogen transmission (Titus et al., 2006; Kazimirova and 
Stibraniova, 2013; Wikel, 2013) Hard ticks have been divided 
into two main phylogenetic lineages, the prostriates (contain- 
ing the genus Ixodes) and the metastriates (containing all other 
hard tick genera). The divergence between these lineages is esti- 
mated to have occurred 124 million years ago (range from 101 to 
166) (Mans et al., 2012). While salivary proteins from prostriate 
and metastriate lineages differ in amino acid composition, they 



accomplish similar goals in terms of inhibiting host responses 
to allow successful blood feeding. However, some differences 
are likely based on differential vectorial capacity. In this study, 
cutaneous host responses to D. andersoni were measured in an 
effort to understand host factors effecting pathogen transmis- 
sion. Affymetrix GeneChip Mouse Genome 430A 2.0 arrays were 
used to measure host gene expression in the skin of D. ander- 
soni infested mice at 12, 48, 96, and 120 h post- infestation (hpi) 
during primary infestations and 120 hpi during secondary infes- 
tations. In addition, histopathological analysis of bite-site lesions 
from primary infestation time points and quantitative real-time 
PCR analysis of lymph nodes from secondary infestation time 
points were analyzed. These studies allow us to describe the cuta- 
neous host response during primary and secondary infestations, 
measure changes in gene expression patterns across time, view 
potential patterns related to tick immuno-suppression, correlate 
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the histopathology analysis to the gene expression data, and 
measure systemic responses in the draining lymph node. These 
analyses are important for understanding the context of pathogen 
transmission and tick rejection at the tick-host inGAterface. 

METHODS 
TICKS 

Pathogen-free D. andersoni colonies were maintained in our lab- 
oratory as described (Bouchard and Wikel, 2005; Heinze et al., 
2012a,b). All life cycle stages were kept in sterile glass vials with 
mesh tops in desiccators at 22°C containing saturated solutions 
of KNO3 to obtain desired relative humidity with a 16:8 h pho- 
toperiod. For routine colony maintenance adult and/or nymphal 
ticks were fed on New Zealand white rabbits and nymphs and/or 
larvae were fed on mice. 

ANIMALS 

BALB/c mice used in this study were obtained from The Jackson 
Laboratory (Bar Harbor, ME). Mice were cared for in accor- 
dance with a protocol approved by the Institutional Animal Care 
and Use Committee (IACUC) of the University of Texas Medical 
Branch. 

TIME COURSE INFESTATIONS 

To perform time course infestations, 6-8 week old female BALB/c 
mice were placed into individual restrainers or anesthetized with 
a 150 u,l intraperitoneal injection containing 10 mg/mL ketamine 
(Fort Dodge Animal Health, Fort Dodge, IA) and 1 mg/mL 
xylazine (Phoenix Pharmaceutical, St. Joseph, MO) in sterile 
PBS (Gibco, Life Technologies, Carlsbad, CA) and infested with 
pathogen-free nymphal ticks. Ticks were allowed to attach for 
approximately 1 h and unattached ticks were discarded. Mice 
were then removed from restraints and housed individually. 
Secondary infestations involved two rounds of infestation. Mice 
were infested with nymphal ticks that were allowed to complete 
their feeding cycle (4-5 days). 14 days after the last primary 
infestation tick completed feeding, mice were re-infested with 
nymphal ticks using the same protocol described above. Bite sites 
were analyzed at 12, 48, 96, and 120 hpi during primary infesta- 
tions, and at 120 hpi during secondary infestations. Three mice 
were measured at each time point in all micro-array or PCR-array 
experiments; controls consisted of 3 similarly housed but tick- free 
mice. Skin biopsies (4 mm) at the tick feeding loci (around ear) 
were harvested and stored in RNAlater (Ambion) at — 20°C for 
RNA extraction or fixed in 10% neutral buffered formalin for 48 h 
and then stored in 70% ethanol at 4°C for histopathology. For 
histopathology, skin biopsies were harvested with ticks attached. 

RNA ISOLATION 

Ticks were removed from all skin biopsies before RNA extraction. 
Tissue samples were homogenized individually in 1 mL Trizol 
(Life Technologies, Carlsbad, CA) using an Ultra-Turrax T8 (Ika, 
Wilmington, NC) tissue disperser. The Trizol protocol was fol- 
lowed until phase separation. The aqueous phase was retained 
and one volume 70% ethanol (Acros Organics) was added and 
the samples applied to RNeasy micro kit (Qiagen) columns. The 
RNeasy protocol was then followed, including the in-column 



DNase digestion step. One extra wash step with buffer RW1 and 
two extra wash steps with buffer RPE facilitated the isolation 
of high quality RNA. All samples were eluted in nuclease-free 
water. After extraction, RNA was quantitated spectrophotomet- 
rically using a NanoDrop ND-1000 (NanoDrop Technologies, 
DE). All samples were required to read greater than 1.8 on 
both A260/A28O and A260/A230 ratios. For subsequent microarray 
analysis, quality of the purified RNA was assessed by visualiza- 
tion of 18 and 28 S RNA bands using an Agilent BioAnalyzer 
2100 (Agilent Technologies, CA). Resulting electropherograms 
were used in the calculation of the 28.18 S ratio and the RNA 
Integrity Number, which was greater than 6.8 in all samples 
(Schroeder et al., 2006). For subsequent real-time PCR analy- 
sis, RNA integrity was determined by denaturing (formaldehyde) 
agarose gel electrophoresis followed by staining with Sybr Gold 
stain (Invitrogen). Visualization of clear ribosomal bands indi- 
cated minimal degradation. Eluted RNA samples were aliquoted 
and stored at — 80°C until used. 

HOST GENE EXPRESSION PROFILING USING AFFYMETRIX GENECHIPS 

Total RNA (500 ng) was converted to cRNA for microarray analy- 
sis using the Ambion MessageAmp™ Premier RNA Amplification 
Kit (Life Technologies Corporation, CA) according to man- 
ufacturer's instructions. Total fragmented cRNA (10 tig) was 
hybridized to the Affymetrix GeneChip Mouse Genome 430A 2.0 
array according to the manufacturer's (Affymetrix, CA) condi- 
tions. The chips were washed and stained in a GeneChip Fluidics 
Station 450 and fluorescence detected with an Affymetrix-7G 
Gene Array scanner using the Affymetrix GeneChip Command 
Console software (AGCC1.1). The resulting CEL files were 
uploaded to iReport™ (Ingenuity Systems), and the default data 
analysis was utilized. The resulting values were then filtered for 
p-values < 0.05 and a fold change < —1.5 or > +1.5 and these 
data are reported as significant findings (Figure 1). 

GENE ONTOLOGY(GO) ANALYSIS 

Broad patterns of significantly modulated genes were visualized 
using heatmaps plotted in the R package gplots (Warnes et al., 
2013). Two gene lists were used to generate these maps; the 
first included all genes modulated during the primary infesta- 
tion and the second included all genes modulated at any time 
point in the study (Figure 2). These heat maps suggested mod- 
ulated genes grouped into early primary infestation (12-96 hpi), 
late primary infestation (120 hpi), and secondary infestation 
patterns (120 hpi). The significantly modulated genes in these 
groups were compiled and used to generate a Venn diagram 
(Figure 3) through an online interactive tool (Oliveros, 2007). 
The lists of unique and shared genes in each group identi- 
fied in the Venn diagram were analyzed using the Database for 
Annotation, Visualization, and Integrated Discovery (DAVID) 
(Huang da et al., 2009a,b) website using the Affymetrix GeneChip 
Mouse Genome 430A 2.0 array as a background. The functional 
annotation clustering tool was used to group similar gene ontol- 
ogy terms into clusters with a mathematically-generated ranking 
score. These clusters can then be named based on the gene ontol- 
ogy terms in each cluster, giving a concise over- view of the results 
(Table 1). 
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FIGURE 1 | Number of up and downregulated genes at each time point 
in the microarray data. Mice were infested with D. andersoni nymphs and 
4 mm skin biopsies collected at 12, 48, 96, and 120 hpi during primary 
infestations and 120 hpi during secondary infestations. Gene expression 
profiling was measured using Affymetrix mouse genome 430A 2.0 arrays. 
Significance was assessed using iReports data analysis methods, and 
results were filtered for genes with fold changes greater than ±1.5 and 
p-values less than 0.05. DAR D. andersoni primary; DAS, D. andersoni 
secondary; hpi, hours post infestation. 



VALIDATION OF ARRAY DATA 

The GeneChip results were validated by an additional experi- 
ment. A separate set of mice were infested with ticks and tissues 
harvested as described above. For the D. andersoni microar- 
ray experiment, gene targets (Table 2) were validated at 12, 48, 
96, and 120 hpi during primary infestations and 120 hpi dur- 
ing secondary infestations using 3 mice/time point. Primers 
were purchased from Integrated DNA Technologies; these primer 
sequences are provided in Table 3. Primers were mixed with RT 2 
SYBR green qPCR master mix (Qiagen) and aliquoted into iCy- 
cler iQ PCR plates (Bio-Rad) using an epMotion 5075 automated 
pipetting system (Eppendorf). Plates were sealed and stored at 
— 20°C until use. For each real-time PCR run, the RT 2 First 
Strand Kit (Qiagen) was used to convert 1 ug total RNA into 
cDNA, which was then loaded onto PCR plates using the epMo- 
tion 5075 automated pipetting system (Eppendorf). In some 
cases, the epMotion system was not functioning properly and 
the plates were loaded with primer and/or cDNA using an 8- 
channel pipette. These plates were run on an iCycler iQ5 real-time 
PCR instrument (Bio-Rad) with the following cycling protocol: 
lOmin at 95°C; 15 s at 95°C, 1 min 60°C for 40 cycles, and an 
80-cycle (+0.5°C/cycle) 55-95°C melt curve. Every run included 
hypoxanthine guanine phosphoribosyl transferase (Hprt) and 
heat shock protein 90 alpha (Hsp90abl) as endogenous control 
genes, and "no template" and "no first strand" controls. HTqPCR, 
an R-based program for real-time PCR data analysis (Dvinge and 
Bertone, 2009), was used to analyze data using the delta-delta 
Ct method for gene expression normalization and measurement, 
and the linear models in microarray analysis (LIMMA) package 
for statistical comparisons between infested and tick-free mice. 




FIGURE 2 | Heatmaps showing changes in gene expression across 
time. (A) Changes in genes significant at any time point in the microarray 
study. (B) Changes in genes significant during the primary infestation only. 



QUANTITATIVE REAL-TIME PCR ANALYSIS OF MOUSE LYMPH NODES 
AFTER D. ANDERSONI FEEDING 

Draining cervical/auricular lymph nodes were harvested from 
mice used for validation of the D. andersoni microarray data set 
during secondary infestations at 120 hpi. Quantitative real-time 
PCR analysis was undertaken using the gene targets in Table 4. 
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FIGURE 3 | Venn diagram showing the overlap of significantly 
modulated genes between time points in the microarray study. 

(A) Venn diagram for primary infestation time points. (B) Venn 
diagram for early time points (DAP12, 48, and 96 hpi), later time 



These results were compared against lymph nodes from tick-free 
mice. RNA was isolated and PCR run as described above. 

RESULTS AND DISCUSSION 

The number of up and downregulated genes meeting our filter- 
ing criteria at each time point is shown in Figure 1. Although 
there was a small decrease in the number of upregulated genes 
between 12 and 48 hpi, the general trend was an increasing num- 
ber of upregulated genes across time. The downregulated genes 
show a different profile, with more genes modulated at early and 
late time points than those in the middle. A list of all genes mod- 
ulated at any time point (Supplement data 1) in the study was 
used generate a heatmap (Figure 2A). This heatmap suggested a 
similar gene expression profile for the first three time points (12, 
48, and 96 hpi during primary infestations) that then shifted to a 
second pattern that was similar between 120 hpi in primary and 
secondary infestations. This suggested there was a shift in gene 
expression profiles between early and late time points, as well as 
between primary and secondary infestations. 

To provide a higher-resolution analysis of the primary infes- 
tation time points, a second heatmap was constructed using 
only the significantly modulated genes in the primary infestation 
(Figure 2B). While this heatmap shows some differences between 
12, 48, and 96 hpi, the overall similarity seen in the first heatmap 
was preserved. In addition, the dramatic change in expression 
profile between the first three time points and the 120 hpi primary 
infestation time point was maintained. 

A Venn diagram for the primary infestation time points 
(Figure 3A) shows more uniquely modulated genes at 12 hpi and 
120 hpi than at other time points, again suggesting differences 
between early and late responses. To explore the patterns observed 
in the heatmaps, all significantly modulated genes from 12, 48, 
and 96 hpi primary infestation were combined into a single data 
set. This list, the list of modulated genes at 120 hpi primary 
infestation, and the list of modulated genes at 120 hpi secondary 
infestation were used to create a Venn diagram (Figure 3B). This 
revealed 137 genes unique to early primary infestations, 131 
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points (DAP120 hpi), and secondary infestation (DAS120 hpi). The Venn 
diagram in B was constructed based on variations in gene expression 
profiles suggested in the heatmaps. All Venn diagrams were created 
using (Oliveros, 2007). 



genes unique to late primary infestations, 651 genes unique to 
late secondary infestations, and 57 genes shared across all time 
points. These lists, along with the lists of genes shared between 
each group were submitted to the Database for Annotation, 
Visualization, and Integrated Discovery (DAVID) bioinformat- 
ics database using the Affymetrix Mouse Genome 430A 2.0 array 
as a background list. The functional annotation-clustering tool 
was used to group similar ontology terms together, and the top 
ten clusters containing significant ontology terms (p < 0.05) were 
named based on the terms included (Table 1). 

THE CORE RESPONSE (GENES SHARED BETWEEN ALL TIME POINTS) 

Analyzing the shared response between all time points allows the 
description of a "core" response to tick feeding. Genes in the 
core response were significant at early, primary, late primary, and 
late secondary. Significant gene ontology clusters were chemo- 
taxis and inflammation, cation homeostasis, intermediate fila- 
ments and keratinization, carbohydrate binding, and regulation 
of oxidatives stress. 

Genes in the chemotaxis and inflammation clusters 
were chemokines, cytokines, and anti-microbial molecules. 
Chemokines Ccrl, Ccl2, Ccl6, Ccl7, Ccll2, Cxcll, Cxcl2, and Pf4 
(platelet factor 4 or Cxcl4) were upregulated and are consistent 
with monocyte and neutrophil migration into the bite site. 
Additional chemoattractants for neutrophils such as the SI 00 
molecules (calgranulins) and Saa3 (serum amyloid A3) were 
also upregulated. Antimicrobial proteins such as Lcn2 (lipocalin 
2) and Ptx3 (pentraxin related gene) that function in iron 
sequestration and opsonization respectively were upregulated. 
These results combined with the upregulation of IL-lb suggest a 
potent innate-like immune response dominated by macrophages 
and neutrophils is a hallmark of the anti-tick response. 

Gene ontology terms in clusters 3 and 6 were related to cations 
or metal ions. Genes in these clusters related to ion regula- 
tion were Mtl and Mt2 (metallothionein); these molecules are 
involved in regulating copper and zinc, but also are involved 
in reactive oxygen species scavenging, cell transcription, and 
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Table 1 | Top ten clusters from DAVID analysis based on lists from 
Figure 3B Venn diagram. 



Core DAP1 2-48-96 only 



1 


Chemotaxis, immune 


1 


Response to DNA damage and 




response, response to 




cellular stress, DNA repair 




wounding 






2 


Chemotaxis, cell migration 


2 


Intracellular 








non-membrane-bounded 








organelle 


3 


Cation homeostasis 


3 


Nucleotide and ATP binding 


4 


Intermediate filament, 


4 


Regulation of transcription, DNA 




Cy lUb K.C It; IUI I U I y d I I IZd LIU 1 1 dl IU 




Ull iu 1 1 iy 




structure 






5 


Keratinization, keratinocyte 


5 


Circadian rhythm 




and epidermal cell 








differentiation 






6 


Response to metal ion 


6 


Chromosome, chromatin binding 








and organization 


7 


Carbohydrate binding 


7 


Nuclear body, nucleoplasm 


8 


Regulation of blood pressure 


8 


mRNA metabolic process 


9 




9 


Regulation of transcription 


10 




10 


Cardiac muscle contraction 


DAP120-DAP12-48-96 


DAP120 only 


1 


No clusters with significant 


1 


Regulation of cell and 




terms 




cell-substrate adhesion 


2 




2 


Blood vessel development 


3 




3 


Embryonic development 


4 




4 


Nuclear or intracellular organelle 








lumen 


5 




5 


RNA processing 


6 




6 


Cell growth 


7 




7 


MAP kinase activity 


8 




8 


Positive regulation of transcription 








and RNA metabolic process 


9 




9 


Meisis (cell cycle) 


DAP120-DAS120 


DAS120 


1 


Extracellular region 


1 


Inflammatory and defense 








response 


2 


Proteinaceous extracellular 


2 


Leukocyte and lymphocyte 




matrix 




activation, T cell differentiation 


3 


Cell surface 


3 


Positive regulation of immune 








response, immune 








response-regulating cell surface 








receptor signaling pathway 


4 


Response to wounding 


4 


Lysosome 


5 


Cellular homeostasis 


5 


T cell activation, leukocyte and 








lymphocyte proliferation 


6 


Fatty acid biosynthetic process 


6 


Chemotaxis 


7 


Regulation of phosphorylation, 


7 


Regulation of cytokine production 




protein kinase activity 






8 




8 


Regulation of immune cells 


9 




9 


Regulation of adaptive immune 








response 


10 




10 


Endocytosis 



(Continued) 



Table 1 | Continued 
DAS120-DAP12-48-96 



1 Keratinocyte, epidermal cell 
differentiation 

2 Rhythmic process 

The core analysis represents significantly modulated genes shared between 
early primary (DAP12, 48, 96), late primary (DAP120), and secondary infestations 
(DAS120). DAP, D. andersoni primary; DAS, D. andersoni secondary. 



Table 2 | Gene list for quantitative real-time PCR validation of the D. 
andersoni microarray study. 



Arg1 


Cxcl5 


il4 


S100a9 


Ccl12 




II6 


sele 


Ccl6 


Hprt 


Irak4 


Serpinel 


Ccl7 


Hsp90ab1 


krt16 


Socs3 


Ccr1 


ifng 


krt6b 


Stat3 


Chi3l1 


il 10 


Map3k6 


Wnk1 


CledOa 


il12a 


Mt1 


No RT 


Clec4e 


111 b 


ptges 


No temp 



immune responses. Studies with metallothionein knockout mice 
suggest Mtl and 2 play protective roles in inflammatory set- 
tings (Inoue et al, 2009; Thirumoorthy et al., 2011). Xanthine 
dehydrogenase, an enzyme involved in the oxidative metabolism 
of purines, was also upregulated. Xanthine dehydrogenase may 
be converted to xanthine oxidase by sulfhydryl oxidation, and 
xanthine oxidase can be an important source of reactive oxygen 
species (Sakuma et al., 2012). Thus the generation and scaveng- 
ing of reactive oxygen species appears to be the primary function 
of the genes in these clusters. 

Clusters 4 and 5 were dominated by gene ontology terms 
related to intermediate filaments and keratinization. Indeed, 
keratin intermediate filaments Krt6a, 6b, and 16 were upregu- 
lated. Krtl6 has been shown to pair with either Krt6a or 6b to 
form essential intermediate filaments in epithelial wound repair. 
Specifically, these intermediate filaments are required to support 
keratinocyte migration into the wound site (Paladini et al., 1996; 
Wong and Coulombe, 2003). In addition to these molecules, 
Sprrlb (small proline rich lb) was upregulated. The encoded 
protein is a significant component of the cornified cell enve- 
lope that maintains a permeability barrier across epithelial tissues 
(Patel et al., 2003). Downregulated genes included a minor colla- 
gen (type 11 alpha 1), the cytoskeletal component spectrin, and 
a little-studied molecule Limchl (LIM and calponin homology 
domains 1). These results suggest a keratin intermediate filament- 
based wound healing response is activated at the tick feeding 
lesion, however it is surprising that other cytoskeletal elements are 
downregulated, a pattern that is even more striking in the early 
primary infestation as discussed below. Thus keratin intermedi- 
ate filaments but not other cytoskeletal molecules are activated in 
response to tick bites. 

Gene ontology cluster 7 was related to carbohydrate bind- 
ing. The primary molecules in this group were c-type lectin 
receptors Clec4d and Clec4n, although other molecules such as 
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Table 3 | List of gene targets and primer sequences used in the 
D. andersoni quantitative real-time PCR validation of the microarray 
data and also in quantitative real-time PCR analysis of lymph nodes 
at 120 hpi during secondary infestations. 



Gene name 


Primer 1 (5'-3') 


Primer 2 (5'-3') 


RORgt 


L^uyuiydydyyyuiLUdu 


Tmnammantammnn a nattan 

i yudyydy Ldyy uudud LLdu 


Gapdh 


yLyydyiudLduiyyddUdLyLdy 


Aatmmtmaammtnmmtmtn 

™ Ly y Lyddy y Luyy Ly Ly 


Tbx21 


ijaayaijijaija lUUabaaoUa LU 


Ttna anna ninannania nam 
i Luaauuay uauuayauay 


Arn1 

/-u y i 


rt y i y LL y ai y LljC, y L y Ly^y^ 


(^aatnnaanantrantntrint 


Stat4 


f a a nt nntntm tn a pr a t nit 
odalj LUU LU iy LUaUUa Ly I 


f m tn a a a rt nt a t nit nn a m tm a 
L>y Luaaay u La Ly Luuay Ly a 


Chi3l1 


f na tn a a a n r* n a t a a m a a c nr 
L^ua LUaaayuua Laayaauyu 


nnamaaanannaanntmaama 
uuay aaauauuaauu Lyaaya 


Hdc 


f"nannmaatnanaaannanam 
ucHjUycia ILaLaaciUbaUay 


tntanntnnmanatmnnaa 
lu Lauu Luuy aua Ly uuaa 


Irak4 


fanpantanttnannttrsm 
L>ayuayiay LLy ay y llu auy 


Anannnaaatntmanatntana 
n\uauuuaaa lu Ly aua lu Laua 


Statl 


ni a nttn a n a c a n a m a a a tn a a nt n 
y do llu ay auauay aaa Luaau lu 


Ttmanaaamannanmnntt 
i Ly auaaayauuauy uu ll 


Stat6 


amttnttnntmnttnnmatm 
dy LLL.LLUL.LytjLLUL.yd Ly 


u nna nnatn amanaaatanttn 
LDUUdUUd LUdydUddd LdU LLU 


Tnfhl 
i yiu i 


nm tn m a m ttt mtt a tnttt n ptn 

uy Lyy dy LLLy LLd lu LLiyu Ly 


(A anmtnantmma mttmtanm 

uduy LUdu iy ydy LLy Lduy 


^tatT 
o LdLo 


Agctcctcag tcacgatca 


gttcaagcacctgacccttag 




f*ammaatmamtananmaamnn 

Ljdyy dd iydy LdUduy ddyuu 


f^t attan am m amtmannnamatm 

ld id LLdUdy ydy lyduuudyd iy 


1112a 


ntntnmttnttntmtanttnna 
uiu Luy llu LLy iy Lay LLUua 


anamatmanatmmtmaamanm 
auay a Ly aua iy y Ly aay auy 


ruxpo 


C tgtcttccaag tctcg tctg 


Ctggtctctgcaggtttagtg 


Lcn2 


f*ntntm natatttnnna m a nit 
LjU Ly Ly Ud Ld LLLUUUdydy L 


f*tanaatmtnanntnnatnntm 
L> LdUdd Ly LUdUU LUUd LUU Ly 


I1 1 h 
II I u 


Cgagatttgaagctggatgc 


Tgacagtgatgagaatgacctg 


1117a 
1 1 I /d 


gagcttcccaga tcacagag 


Agactacctcaaccgttcca 


1110 
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Ccl7, Pf4, and Ptx3 were also present. C-type lectins are expressed 
on dendritic cells and macrophages, and play important roles 
in orchestrating the inflammatory response. Their importance 
host responses to tick feeding is strongly suggested by studies 
showing the I. scapularis salivary protein Salp 1 5 can bind and acti- 
vate DC-sign (CD209) (Hovius et al, 2008). Clec4n (dectin-2) is 



Table 4 | Gene list for quantitative real-time PCR analysis of mouse 
lymph nodes after secondary infestation with D. andersoni for 120 h. 
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thought the associate with FcRy and signal through a Syk-Card9 
pathway to activate NFkB (Geijtenbeek and Gringhuis, 2009). 
Ligands include mannose, fungal antigens, and dust mites (Patel 
et al., 2003), although others likely exist. The precise role of c- 
type lectins in tick feeding is unknown, but these results suggest 
they may act as important pattern recognition receptors and may 
modulate downstream responses. 

The final cluster contained terms related to regulation of blood 
pressure. This primary gene in this group, Wnkl (lysine deficient 
protein kinase 1), is downregulated and it influences blood pres- 
sure by regulating salt absorption in the kidney. Wnkl also plays 
a similar role in sweat glands (Kahle et al, 2008). Upregulated 
genes in this group include Hmoxl (heme oxygenase 1), and 
Gchl (GTP cyclohydrolase 1). Hmoxl can play a role in pro- 
tecting tissues from oxidative damage, while Gchl can produce 
either NO or reactive oxygen species depending cofactor availabil- 
ity. Thus the gene ontology terminology may be misleading here, 
as the genes appear more related to regulating oxidative status 
than blood pressure. 

EARLY PRIMARY INFESTATION 

Gene ontology analysis of the early primary infestation revealed a 
number of clusters related to nucleotide processing such as DNA 
repair, nucleotide binding, transcription, and mRNA metabolism. 
A review of the genes represented by these clusters suggested 
they could be broadly grouped into DNA repair molecules, DNA 
helicases of the SWI/SNF family, transcription factors, compo- 
nents of the spliceosome, and mRNA metabolism. DNA helicases 
of the SWI/SNF family are not active helicases, but function 
as DNA translocases thought to modulate chromatin structure 
(Durr et al., 2006). Thus most of these molecules can be related 
to transcription, from chromatin remodeling, transcription fac- 
tor binding, RNA splicing, and mRNA metabolism. Interestingly, 
nearly all of these genes were downregulated, suggesting that tran- 
scription is decreased early in the primary infestation. Of the 
few upregulated transcriptional regulators, Nfkbia and Tsc22d3 
were of particular interest because they have been shown to 
inhibit NFkB and AP-1 pro-inflammatory pathways (Beaulieu 
and Morand, 2011). DNA repair molecules were also downregu- 
lated. In contrast, molecules in the circadian rhythm cluster, such 
as Perl were upregulated. Overexpression of Perl in tumor cells 
was shown to increase DNA-damage induced apoptosis and sig- 
nificantly decreased cell growth in untreated cells (Gery et al., 
2006). This suggests some connection between DNA repair, cir- 
cadian rhythm molecules, and growth arrest in the skin at the 
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bite site may exist. Terms in the final cluster were related to 
cardiac muscle contraction, although review of these genes and 
genes in the highly significant (though un-clustered) GO term 
actin binding (GO:0003779) revealed molecules related to mus- 
cle contraction, actin binding, and cytoskeletal adaptor molecules 
that may relate signaling events to changes in the cytoskeleton. 
Cytoskeletal elements have been shown to play vital roles in 
epithelial wound healing processes by aiding cell migration and 
wound contracture (Abreu-Blanco et al., 2012). Thus the early 
host responses were characterized by a reduction in transcription, 
DNA repair, inhibition of inflammation, potential inhibition of 
the cell cycle, and downregulation of cytoskeletal elements. These 
effects suggest that tick feeding may inhibit early aspects of the 
wound healing response. 

When the genes shared between early and late primary infes- 
tation were submitted to DAVID, no significant clusters were 
observed. 

LATE PRIMARY INFESTATION 

The first three clusters (regulation of adhesion, blood vessel 
development, embryonic development) and clusters 6 and 7 
(cell growth, MAPK activity) share similar genes and are treated 
together here. In the process of normal epidermal wound healing, 
repair must take place in the epidermis, dermis, and vasculature to 
restore normal function. Particularly for the re-epithelialization 
of the wound, keratinocytes must loose some of their adhesive- 
ness and migrate into the wound area. This process shares similar 
characteristics to the epithelial-mesenchymal transition (EMT) 
(Nakamura and Tokura, 2011) during embryonic development. 
In support of a similar process at the tick bite site, genes known 
to be involved in EMT such as Tdgfl (teratocarcinoma-derived 
growth factor 1 or Cripto-1), Cyr61 (cysteine rich protein 61 or 
CCN1), and Smad5 (SMAD family member 5, a TGF-p" signal- 
ing intermediate) were upregulated. While not specific to EMT, 
the MAPK signaling pathway is known to be important (Leopold 
et al, 2012), increasing the relevance of cluster 7. Additional 
upregulated genes that may be involved were Tnfrsfl2s (tumor 
necrosis factor receptor superfamily, member 12a), Junb (Jun-B 
oncogene), and Epgn (epithelial mitogen). While a specific role 
for the encoded proteins in EMT has not been delineated, these 
molecules may be related to TNF-ot, AP-1, and growth factor 
responses that are known to be involved (Leopold et al., 2012). 

Downregulated molecules included extracellular matrix com- 
ponents Collal (collagen type I alpha 1) and Lamb2 (laminin 
beta 2), calcium ion channel or signaling molecules Tesc (tescal- 
cin) and Pkd2 (polycyctic kidney disease 2), integrin ligand Edil3 
(EGF-like repeats and discoiding I-like domains 3), and growth 
factor related molecules Pdgfrb (platelet derived growth factor 
receptor beta) and Tgfb3 (transforming growth factor beta 3). 
In the final analysis, a convincing case for EMT in wound heal- 
ing responses during late primary infestations cannot be made. 
However, changes in adhesion, cell growth, MAPK activity, and 
the appearance of genes normally associated with embryonic 
development such as Tdgfl are suggestive of a related process. 

The remaining clusters were nuclear or organelle lumen, RNA 
processing, positive regulation of transcription, and meiosis or 
cell cycle related. The majority of the genes in these clusters 



code for RNA binding proteins, RNA processing proteins, and 
transcriptional activators/transcription factors. An interesting 
theme among the RNA-related molecules was an interaction with 
rRNA, implying some interaction with ribosomes. Most of the 
genes in these clusters were upregulated, contrasting with the 
transcription-related genes early in the primary infestation. Thus, 
the data suggests transcription is inhibited during early host 
responses, but then activated near the end of the feeding cycle. 

Genes shared between 120 hpi primary and 120 hpi secondary 
infestation clustered into a number of broad categories such as 
extracellular region, cell surface, cellular homeostasis, and regu- 
lation of phosphorylation. Two clusters were more specific, how- 
ever. Fatty acid biosynthesis contained three genes, Fcerla (IgE 
Fc receptor), Scd3 (stearoyl-coenzyme A desaturase), and Ch25h 
(cholesterol 25-hydroxylase). Fcerla and Ch25h were upregu- 
lated, while Scd3 was downregulated. Fcerla encodes the alpha 
chain of the FceR, a receptor that binds IgE with high affin- 
ity. This receptor is activated by antigen-induced cross-linking, 
and induces many downstream responses including the release 
of histamine and the production of prostaglandins D2 and E2 
(Ugajin et al., 2011). Supporting this role, genes encoding pro- 
teins involved in prostaglandin E synthesis are upregulated at 
120 hpi primary infestation, although they are not shared with 
120 hpi secondary infestation. 

Ch25h encodes a protein that catalyzes the formation of 25- 
hydroxycholesterol (Lund et al., 1998), a molecule that can act as a 
chemoattractant for migrating cells but also inhibits the function 
of B cells (Diczfalusy et al, 2009) and can promote the survival 
of Listeria monocytogenes infected cells (Zou et al, 2011). Scd3 
is an important enzyme in the synthesis of monounsaturated 
fatty acids, many of which are subsequently used as components 
of membrane phospholipids. Interestingly, decreasing fatty acid 
synthesis inhibited rift valley fever virus replication in mouse 
embryonic fibroblasts (Moser et al., 2012), suggesting fatty acid 
synthesis may aid viral replication. These genes suggest a role for 
fatty acid metabolites and the control of fatty acid synthesis at the 
tick bite site. Finally, the response to wounding cluster contained 
pro-inflammatory genes such as 11-6, Cxcl5, and P2ryl2 (puriner- 
gic receptor P2Y) and extracellular matrix molecules Cd44 and 
Timp3 (tissue inhibitor of metalloproteinase 3). These genes are 
involved in inflammation and cell migration, although they are 
often associated with acute inflammation rather than appearing 
late in the induction of the immune response as they do here. 

In summary, late in the primary infestation the innate immune 
response remains the major player. Keratin-based wound healing 
responses and the suggestion of an EMT-type transition driving 
wound re-epithelialization are present. In addition, transcription 
appears to be potentiated in contrast to earlier in the primary 
infestation. Genes shared with the secondary infestation sug- 
gest fatty acid metabolism and the response to wounding are 
important in host responses near the end of the feeding cycle. 

SECONDARY INFESTATION 

There were 651 genes modulated only in the secondary infesta- 
tion. Gene ontology analysis of this group of genes resulted in 
many clusters related to the immune response. These clusters 
shared many genes with each other and they will be discussed in 
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relation to what appeared to be the most distinctive feature(s) of 
each cluster or group of clusters. 

The first cluster contained GO terms related to inflamma- 
tory and defense responses. Genes in this group supported four 
primary processes: complement (synthesis), coagulation (path- 
way synthesis), Toll-like receptor/IL- 1 response, and acute phase 
response. Clusters 2, 5, 8, and 9 were leukocyte and lympho- 
cyte activation, T cell activation and leukocyte proliferation, 
regulation of immune cells, and regulation of adaptive immune 
response, respectively. These clusters shared similar genes, and 
supported the importance of protein kinases (SRC family, Syk, 
and Jak), T-cell costimulation (CD40, CD86, CD3d, and Lck), 
cytokines, immunoreceptors (cytokine, antibody, C-type lectin, 
and Toll-like receptors 1 and 4), and signaling intermediates 
in immunoreceptor signaling. Cluster 3, positive regulation of 
immune response, contained many similar features. However, two 
genes strongly related to iron handling in the innate immune 
response, Slcllal (solute carrier family 11 member 1) and Hpx 
(hemopexin) were upregulated. Cluster 4 contained GO terms 
related to lysosomes. The largest group of upregulated molecules 
in this cluster were genes encoding many different cathepsins. 
Cathepsins are lysosomal proteinases that function primarily in 
protein turnover in the lysosome (Kuester et al, 2008). However, 
additional roles for cathepsins include (but are not limited to) 
antigen presentation (Yamamoto et al., 2012), keratinocyte dif- 
ferentiation (Yamamoto et al., 2012), itch responses (Andoh et al., 
2012), and wound healing (Akira and Takeda, 2004). The chemo- 
taxis cluster contained chemokines Ccl24, Ccll7, Ccl5, Ccl4, Ccl8, 
Ccl3, and Cxcll6 consistent with the migration of lymphocytes, 
monocytes, neutrophils, eosinophils, and basophils into the bite 
site. Additional molecules that supported cell migration such 
as integrins, chemoattractants, and intracellular molecules that 
regulate changes in cell shape were also upregulated. Cluster 7 
contained terms related to regulation of cytokine production. 
These genes code for proteins that function as immunorecep- 
tors, cytokine signaling intermediates, and cytokines. Finally, the 
endocytosis cluster contained genes encoding classical endocy- 
totic proteins such as Clta (clathrin, light polypeptide) and picalm 
(phosphatidylinositol binding clathrin assembly protein). Other 
genes in this cluster were antibody Fc receptors Fcerlg, Fcgr3, and 
Fcgr2b and lectin recepors Clec7a, CD209a CD209d, CD209e, 
and Mrcl. These receptors may function in endocytosis by bind- 
ing antibodies or carbohydrate moieties and triggering endocy- 
totic machinery. This may be an important method of delivering 
tick-derived antigens to the endosomal/lysosomal compartment 
where they could be processed for antigen presentation. 

These results suggest the importance of complement, coag- 
ulation, innate immune responses, and adaptive immunity in 
the host response to D. andersoni nymphs. In particular, the 
importance of Toll-like receptors 1 and 4 and lectin pattern recog- 
nition receptors appear to be important aspects of the innate 
immune response during secondary infestations. Both TLRs and 
C-type lectins have been shown to be important pattern recog- 
nition receptors for sensing pathogens and initiating and shaping 
innate and adaptive immune responses (Akira and Takeda, 2004; 
Geijtenbeek and Gringhuis, 2009). Both TLR2 and DC-SIGN 
(CD209e) have already been shown to be important in anti-tick 



responses by modulating cytokine responses (Hovius et al., 2008; 
Oliveira et al, 2010). These studies support our results and also 
suggest the importance of these receptors is likely broader than 
previously reported. Within the adaptive immune response, these 
results suggest the capture, processing, and presentation of anti- 
gens is an important process during secondary infestations. In 
particular, cathepsins may be important proteins in this pro- 
cess with possible contributions to wound healing and/or itch 
responses. The array results also support the co-stimulation and 
activation of lymphocytes. 

VALIDATION 

A list of genes was chosen to validate the microarray data 
(Table 2) using quantitative real-time PCR. Genes were chosen 
based on significant fold change in the microarray study that were 
involved in important pathways identified in the gene ontology 
analysis. Trends of up- and downregulation were well preserved 
between the microarray and validation study (Figure 4). In gen- 
eral, more genes had significant fold changes in the microarray 
that the validation for 12 and 48 hpi, and the validation had more 
significant fold changes from 96 hpi onwards. At the 120 hpi sec- 
ondary infestation time point, all the gene targets validated except 
Chi3ll and Wnkl. It should be noted that some genes significant 
only at later time points in the array study became significant ear- 
lier in the validation study (i.e., Argl, Cxcl5), while some genes 
significant early in the array study did not become significant 
until later time points in the validation (i.e., Ccll2, Ccl7, Krt6b, 
Krtl6, and S100a9). Interestingly, the trend of more significant 
gene modulation at 12 than 48 hpi seen in the array data was 
upheld in the validation study. 

In summary, the validation study upheld the microarray study 
with a few minor differences. Samples used for the validation 
study were derived from a completely separate infestation exper- 
iment than the array study, making it an extremely rigorous 
test and likely contributing to the small variations seen between 
data sets. 

LYMPH NODE STUDY 

The host response to ticks has frequently been reported to be 
polarized toward the production of Th2 cytokines. However, in 
our earlier work with I. scapularis, the gene expression profile in 
the skin was consistent with a more mixed response (Heinze et al., 
2012a,b). Consistent with these results, gene expression profiling 
in the skin at D. andersoni bites did not yield a conclusive pic- 
ture of CD4+ T cell polarization. While IL-4 was present, so were 
IL-6 and IL-lb, cytokines more associated with Thl7 responses 
than Th2. To assay the systemic host response, the expression 
of a panel of genes (Table 4) related to CD4+ T cell differen- 
tiation was measured in lymph node samples from naive and 
120 hpi secondary infestation mice. Important transcription fac- 
tors for T regulatory cells (Foxp3), Thl7 cells (Rorgt), Thl cells 
(Tbx21) and Thl cytokine (IL-12a) were either downregulated 
or unchanged between infested and control mice. In contrast, 
the Th2 cytokine IL-4 was upregulated (Figure 5). Thus, the 
results of this experiment strongly suggest the systemic response 
to repeated D. andersoni nymphal feeding in mice is polarized 
toward a Th2 profile. 



Frontiers in Microbiology | Microbial Immunology 



May 2014 | Volume 5 | Article 198 | 8 



Heinze et al. 



Host immune responses to tick feeding 



DAP12hpi DAP48hpi DAP96hpi DAP120hpi DAS120hpi 




Fold change over tick-free mice 

FIGURE 4 | Validation of microarray data by quantitative real-time PCR. LIMMA (linear models in microarray analysis) implemented in HTqPCR, and 
A list of genes significantly modulated in the microarray study (Table 2) R based program for qrt-PCR analysis (see methods chapter). All the 
were validated using quantitative real-time PCR on bite site skin samples genes marked with an asterisk were significant (p < 0.05) as compared to 
from a separate infestation experiment. Significance was assessed using tick-free mice. 



HISTOLOGY 

The bite site is characterized by a large quantity of attachment 
cement secreted by the tick to anchor it to the host (Figure 6). 
The amount of cement increases throughout the feeding period, 
and it contains readily evident layers, suggesting that additional 
cement is secreted in waves throughout the feeding cycle. The 
deepest penetration of the hypostome into the skin is very shal- 
low, in most cases barely passing through the epidermis. At 12 hpi, 
a mild inflammatory infiltrate was present in the dermis, espe- 
cially concentrated at the junction between dermal connective 
tissue and underlying adipose tissue. A small lesion was located 
immediately under the hypostome and likely represents the feed- 
ing "pool." At 48 hpi, the number of inflammatory cells had not 
increased from 12 hpi, mirroring the reduction in gene expression 
seen at this time point. The feeding lesion is very well defined, and 
extravasated erythrocytes are readily evident around the hypos- 
tome. Some thickening of the epidermis is evident. By 96 hpi, 
the inflammatory infiltrate has increased dramatically and the 
feeding lesion appears to have moved deeper into the dermis. 
Epidermal thickening and extravasated erythrocytes and leuko- 
cytes are present. At 120 hpi, most of the changes at 96 hpi are 
intensified. The infiltrate is very dense, the epidermis is markedly 
thickened, the feeding lesion is poorly defined and the dermal 
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FIGURE 5 | Lymph node gene expression profile supports a systemic 
Th2 response. Gene expression changes between secondary infestation 
120 hpi draining lymph nodes and lymph nodes from tick-free mice 
measured by quantitative real-time PCR. Significance was assessed using 
LIMMA (linear models in microarray analysis) implemented in HTqPCR, and 
R based program for qrt-PCR analysis (see methods chapter). All the genes 
shown were significant (p < 0.05). 
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FIGURE 6 | Histological analysis of D. andersoni bites during primary 
infestations. Biopsies from tick feeding sites were fixed in zinc fixative and 
embedded in paraffin. Paraffin blocks were carefully sectioned and slides 
showing the entry of the hypostome into the skin were stained with H&E. In 



the top two rows, relevant structures such as the tick, tick cement, 
epidermis, palps, and feeding lesion have been labeled. The second row 
appears different that the others because the orientation of the tick in the 
block was different. 



tissue near the hypostome appears to be loosing its normal 
architecture. Extravasated erythrocytes and leukocytes are promi- 
nent. In the microarray, the chemotaxis of immune cells into 
the bite site, and aspects of wound healing responses were espe- 
cially prominent. The histological analysis supports these basic 
conclusions, showing an inflammatory infiltrate that increases 
across time, and significant changes in the epidermal and dermal 
compartments near the feeding tick. 

CONCLUSIONS 

The microarray data suggests the chemotaxis of neutrophils, 
monocytes, and other cell types, the production and scavenging of 
reactive oxygen species, keratin-based wound healing responses, 
and C-type lectins are important host responses to tick feed- 
ing conserved across primary and secondary infestations. Early 
time points in the primary infestation imply the inhibition of 



transcription and RNA processing. This inhibitory profile may 
represent tick-induced modulation of early host responses. Later 
time points in the primary infestation provide additional sup- 
port for wound healing activities including the suggestion of 
EMT. The secondary infestation is characterized by increases in 
these same pathways, along with the addition of complement, 
coagulation, Toll-like receptors, acute phase response, antigen 
presentation, and the activation of lymphocytes. These results 
were supported by quantitative real-time PCR validation. The 
systemic murine response to D. andersoni nymphal feeding was 
polarized toward a Th2, based on real-time PCR analysis of drain- 
ing lymph nodes from secondary infestations. The histological 
analysis supported the chemotaxis of neutrophils and monocytes 
into the bite site and the marked thickening of the epithelial layers 
may be the result of the upregulation of genes related to kera- 
tinization and wound healing. In conclusion, during D. andersoni 
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feeding infiltration of inflammatory cells increases across time 
concurrent with significant changes in the epidermal and dermal 
compartments near the feeding tick. The importance of changes 
in the epidermal layer in the host response to ticks is not known, 
however, it is possible the host attempts to "slough off" the tick 
by greatly increasing epithelial cell replication. 
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